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Abstract 

Doubly periodic (periodic both in time and in space) solutions for the 
Lagrange-Euler equation of the (1 + l)-dimensional scalar <y2 4 theory are 
considered. The nonlinear term is assumed to be small, and the Poincare- 
Lindstedt method is used to find asymptotic solutions in the standing wave 
form. The principal resonance problem, which arises for zero mass, is solved 
if the leading-order term is taken in the form of a Jacobi elliptic function. It 
have been proved that the choice of elliptic cosine with fixed value of module k 
(k « 0.451075598811) as the leading-order term puts the principal resonance 
to zero and allows us to construct (with accuracy to third order of small pa- 
rameter) the asymptotic solution in the standing wave form. To obtain this 
leading-order term the computer algebra system REDUCE have been used. 
We have appended the REDUCE program to this paper. 



1 Introduction 



1.1 Periodic solutions of nonlinear equations 

Recently, periodic solutions of the nonlinear wave equation: 

^-^ )= °' (1) 

where denotes a continuous function on 1R such that g(0) = 0, are discussed inten- 
sively. In the 1970s and 1980s existence theorems for a wide class of continues functions 
g{ip) satisfying definite boundary conditions in x have been proved. 
J. M. Coron M has proved that the following problem: 

f d 2 <p{x,t) d 2 ip(x,t) 



dx 2 
lim <p(x,t) 



dt 2 
0. 



-g(<p) = o, 



x,teJR, g(0) = 0, 
x,t eJR, 



(2) 



x,t,T g IR, 



4 ip(x,t + T) = (p(x,t), 

has a T-periodic (in time) solutions only if g'(0) > C^f-) ■ 

In particular, if g(ip) = sm(<p), there are no T-periodic solutions when T < 2tt (for 
any period T > 2-7T an explicit solution is known, see e.g. G. Lamb @). If g((p) = — sin(^) 
or g(tp) = ±y? 3 , then there is no nontrivial solution of problem (2) (for any period T). To 
obtain periodic solutions of massless tp 4 theory one has to change the boundary conditions. 

Generalizing the Rabinowitz theorem ||, H. Brezis, J. M. Coron and L. Nirenberg M 
have proved that for each T, which is a rational multiple of ir, the following problem: 

f #V0M) d 2 ip(x,t) 



dx 2 dt 2 

(p(0,t)=<p(TT,t) = 0, 



g(<p) = 0, x,teJR, g(Q) = 0, 
t e 1R, 
x,t,T e IR, 



[2') 



t <p(x,t + T) = 

has a nontrivial weak T-periodic solution if g(ip) is a continuous nondecreasing function, 
satisfying the following conditions: 

9(f) 



lim 



■DC' 



-jf-gif) 



\<p\-+OC <y2 

3 a, /? > such that VV G IR 
By a weak solution we mean such function tp that 

J J (ip{y u - v xx ) + g(<p)v) dt dx 





g^)diP>(3\g( V )\-a 







for all v G C 2 ([0, 7r] x IR) satisfying the boundary and periodicity conditions. A simplified 
proof of the Rabinowitz theorem for the case 

9(f) = \f\ n -f, 

where n > 0, was given in ||. 
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1.2 Two kinds of doubly periodic solutions 



Our investigation is dedicated to the construction of doubly periodic classical fields in the 
(1 + 1) -dimensional (p A theory. We study the model of an isolated real scalar field (p(x, t), 
described by the Lagrangian density: 

C(<p) = \ t) - <p%{x, t) - MV (:r, t) - E -if\x, t)) . 

The corresponding Lagrange-Euler equation is: 

d 2 (p(x,t) d 2 (p(x,t) , , 9 . o, . 

g x2 ~ g t 2 ~ M M%, t) - e V \x, t) = 0. (3) 

There are two classes of doubly periodic solutions for this equation. If we seek fields 
in the traveling wave form: 

(p(x, t) = ip(x — vt), 

where v is the velocity of the wave motion, then equation (3) reduces to the Duffing's 
equation ||. As is known [0, periodic solutions for this equation are Jacobi elliptic 
functions. Such solutions are well known for both (1 + l)-dimensional || and (3 + ^-di- 
mensional || variants of ip 4 theory. They are used in the Yang-Mills theory |T0[ . 

The second class of doubly periodic solutions consists of functions in the standing 
wave form: 

oo oo 

ip(x, t) = J2Jl C nj sin(n(x - x )) sin(j u (t - *„)), (4) 

n=l j=l 

where x and t are constants determined by boundary and initial conditions. Equation 
(3) is a translation- invariant, so we can restrict our consideration to the case of zero x and 
t without loss of generality. The choice x = corresponds to the boundary conditions 
of problem (2'). The period of f(x,t) in x is taken to be 2n and we are interesting in 
finding the period in t. 

Exact standing wave solutions are not known. Approximate solutions for equation (3) 
with e — 1 and M = 1 have been found under the assumption that all highest harmonics 
are zeros: Vn, j > N : C n j = 0, where N is a large number. Numerical calculations have 
been made [O] for several values of N. The obtained solutions are very close to each 
other, but convergence of the sequence of these solutions as N tends to infinity cannot be 
proved in numerical framework. 



2 Asymptotic solutions in the standing wave form 



Let us consider periodic solutions for equation (3), provided that e is small or, equivalently, 
that these solutions have small amplitudes. 

An asymptotic expansion, containing only bounded functions, is called a uniform 
expansion. The possibility to obtain a uniform expansion, using standard asymptotic 
methods, for example, Poincare-Lindstedt |12| or Krylov-Bogoliubov methods, de- 
pends on values of the frequencies in the zero approximation (for a review of such methods, 
see UlI)- If £ — 0, then equation (3) is linear and has periodic solutions of the form (4) 
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with frequencies in time Qj = \/j' 2 + M 2 , where j G IN. There are two fundamentally 
different cases. 



If 

Vi ,je IN ^ k,ne IN such that 



then it is a nonresonance case and a uniform expansion can be found with any desired 
accuracy using the standard asymptotic methods, for instance, the Poincare-Lindstedt 
method. All secular terms can be eliminated with the proper choice of frequency. 

The resonance case, when there exist two frequencies flj and whose ratio is ratio- 
nal, is more difficult. The Krylov-Bogoliubov method and the standard variant of the 
Poincare-Lindstedt method can be used to find periodic solutions only to a few leading 
orders in e. 

The important example of resonance case is the massless ip A theory when Qj = j 
and all relations of frequencies are rational, which is the principal resonance case. Using 
standard asymptotic methods, one cannot construct a periodic solution even to the first 
order in e. It is possible to transform the differential equations to a system of nonlinear 
algebraic equations in Fourier coefficients and frequencies using the Poincare-Dulac's 
normal form method. The algorithm of this procedure was constructed []I6| , [17]] and 
realized in the software for symbolic and algebraic computation REDUCE [18-20] (the 
system NORT [21, ^]). But an algorithm to solve the resulting algebraic system has yet 
to be created. 



3 The standing wave solutions for the massless ip 
theory 

3.1 Construction of asymptotic solutions via the Poincare- 
Lindstedt method 

The purpose of this article is the construction of standing wave solutions of equation (3) 
with M — 0, using asymptotic methods when £ < 1. We use the Poincare-Lindstedt 
method: introduce the new time t = ut and look for a doubly periodic function <p(x,t) 
and a frequency (in time) u in the form of power series in e: 

<p(x,i,e) 
u(e) 

Expanding the Lagrange-Euler equation in a power series in e, we obtain a series of 
equations. The first two equations of this series are: 

• to zero order in e the equation for (p is: 



oo 
n=0 

oo 

= i + E 

n=l 
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• to first order in e the equation in ip , ip 1 and is: 



^ =2 " 1 ^^ + ^' t) - (6) 

Equation (5) has many periodic solutions. If we select ip (x,t) = sin(x) sin(t), then 
the second equation has no periodic solution, because the frequency of the external force 
sin(3a;) sin(3t) is equal to the frequency of its own oscillations and it is impossible to put 
this resonance harmonic to zero by selecting only uo l . In the massless case to obtain a 
uniform expansion, we must find another doubly periodic solution to equation (5). The 
correct selection of not only the frequency w(e), but also the function tp (x,t), allows us 
to construct the next to leading order term in standing wave form. 



3.2 The condition of existence of periodic solutions 

The general solution, in the standing wave form (4), for equation (5) is the function 

oo 

(p (x, t) = E a n sin(na;) sin(nt) 

n=l 

with arbitrary a n . We have to find coefficients a n so that the function (pi(x, t) is a periodic 
solution for equation (6). If we select (pi(x, t) as a double sum: 



<Pi(x, = E E Kj sin(nx) sin(jt) 

n=l j=l 

with arbitrary b n j, then equation (6) can be presented in the form of Fourier series: 

<9V(M) #ViOM) d 2 ip (x,t) „ ~ N ™™ . 

q x2 ' ~ ' ~ &>i Q i2 ' - t) = E E R nM> b) smH sm(jt) = 

and is equivalent to the following infinite system of the algebraic equations in Fourier 
coefficients of the functions ip (x,t) and (pi(x, t): 



Vra, j : R nj (a.,b) = 0. 



(6') 



This system has a subsystem of the equations in the Fourier coefficients of the function 

ip (x,t): 

V? G IN : %(a) = 0, (7) 

where 



R j:j = 9a j + Saja 3j + 



\ 



oo 2j-l 
6 E ( 2a2 s + a s a 2j+ s) + 3 E a s a 2j-s ~ 32 J 2 ^1 

8 = 1 8 = 1 , 



+ 



+ 3EE a s a p a j+s+p + 3 E E a s a p a s+ P -j 



3-2 3-2 



EE" a ° ■ 



s=l p=l 
s +3 p¥=3 



s=l p=l 
s^j P+j 

P¥=2j-s 



s=l p=l 



s P js-p 
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If the system (7) is solved and all Fourier coefficients of the function ip (x,t) are 
known, then the system (6') transforms into a system of linear equations: R n j(h) = 0, 
which always has only one solution. Hence, we have obtained the necessary and sufficient 
condition of the existence of periodic solutions for equation (6): there exists a periodic 
function <fi(x, t), satisfying equation (6), if and only if the Fourier coefficients of the 
function ip (x,i) satisfy system (7). 

The coefficient a 1 is a parameter determining the oscillation amplitude. In fact, let 
and uj 1 = c^al, then all polynomials -Rjj(a) are proportional to a\: -Rjj(a) = 
a\Rjj(c) and, therefore, the coefficient a x can be selected arbitrarily. 



3.3 The approximate solution of system (7) 

System (7) is very difficult to solve. On the one hand, all Rjj are infinite series and the 
number of equations is infinite too. On the other hand, each equation of this system is 
nonlinear. We restrict ourselves to finding a particular solution. To simplify calculations, 
we assume that the function (p (x,t) contains only odd harmonics. Our goal is to find a 
real solution so we seek c- G 1R. 

To find an approximate solution we apply the Galerkin method: cut off higher diagonal 
harmonics Vj > iV: c 2 j_i = and seek an approximation for if (x, t) in the following form 

(p (x, t) = a x c 2j _ x sin((2j - l)x) sin((2j - l)t) 

Neglecting the higher harmonics of the function (p Q (x,t), we obtain a finite system of 
nonlinear equations, with the number of equations three times as many as that of the 
variables. We solve the system of the N leading equations and substitute the obtained 
values in the other equations. In principle, any system of N nonlinear equations in N 
unknowns can be solved with a computer, using Buchberger algorithm [ 23 , |2~3j [. This 



algorithm has been realized in the standard procedures of the computer algebra system 
REDUCE: the SOLVE operator and the Groebner basis package. This algorithm allows 
us to diagonalize this system of nonlinear equations, constructing an equivalent system, 
which consists of equation in only one variable, equation in two variables and so on. So, 
we can obtain solution of the system, solving only equations in one variable and substitute 
the result into the remaining equations. 

The operating memory and other parameters of real computers are not enough to 
solve very difficult systems. For example, using computer with 128 Mbytes operating 
memory (RAM), we have solved the system, consisting of the iV leading equations of the 
system (7) only for N < 6. On the other hand, even for these values of N the solution of 
system can be found only with finite accuracy (the exact solution can be irrational). 

The SOLVE operator in REDUCE (in "on rounded" mode) rounds off numbers with 
accuracy 5 = 10 -11 . We assume that the system of the N leading equations is solved if 
for all j < N the inequalities \Rjj(c)\ < 5 are true. We also admit that the value of N is 
sufficient to solve system (7) if for all j G IN the inequalities \Rjj(c)\ < 5 are true. 

To find a particular solution of the N leading equations for any N < 50 we have 
written a short program in REDUCE (see Appendix) and obtained that the minimal 
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sufficient value of N for 5 = 10 11 is N = 8 and that the frequency correction is 

u ± = 0.28268003454 a?. 

We also have found numerical values of the fifteen leading Fourier coefficients of 
(f (x,t). The values of Cj are very close to the values of the corresponding terms of 
the following finite sequence: 

d = i d 2j-i = wh ere f 2j _ x = - + q2j ^v = °5 3< 23 )> 

with q = 0.0142142623201. It is easy to verify that substitution of this finite sequence 
gives V j G IN : \Rjj(d)\ < 10~ 12 . 

In other words, the finite sequence d is an approximate solution of system (7). These 
numerical calculations help to find the analytical form of (p (x,t). The following table 
illustrates this interesting result: 



j 


c j 


Ujj {C) 


a j 


n jj\ a ) 


1 

1 


i 
i 


n 

u 


i 
i 


n 

u 


Q 
O 


1 441R9RR1711 v 1 fl~ 2 


O.O A 1U 


1 441R9RR1711 v 1 D~ 2 

l.^lUiUUl 1 11 A 1U 


8 n v 1 n -14 

— O.U A IU 


K 
O 


9 04Q1 71 774D8 v 1 fl -4 

Z.U'ijl i 1 l l ^UO A 1U 


9 1 v 1 n~ 12 

1 A 1U 


9 fl4Q1 71 7741 Q v 1 fl~ 4 


— ^ 4 v i n -15 

— O.^t A IU 


7 


9 Q1974R4Q794 v 1 n~ 6 


7 s v 1 n~ 12 

( .O A 1U 


9 Q1 974fi^1 ^4^ v 1 D~ 6 


Q 7 v 1 n~ 17 
— y. i a iu 


Q 


zL 1 409^41 81 1 ^ v 1 fl~ 8 


O.O A 1U 


4 1409^4^049^ v 1 D~ 8 

^t. L^kVJ^O^kOVJ^k^O A IU 


—9 3 v 1 n~ 18 

— Z.O A IU 


1 1 
1 i 


^ 88^0(^(39014 v 1 n~ 10 


1 n v 1 n~ 14 

1 .U A 1U 


^ 88^0(^07^98 v 1 n~ 10 

O.OOtJUUUU i JZO A IU 


4 Q v 1 n~ 20 

— *±.c/ A IU 


1 3 


8 3fi4881Q9f)7Q x 1 fT 12 


a a x 1 rr 13 


8 36^1 87906^ x 1fT 12 


— Q 7 x 1 rr 22 


15 


1.18901919266 x 10" 13 


-7.8 x 10~ 22 


1.1890496659 x 10~ 13 


-1.8 x 10~ 23 


17 





4.4 x 10" 12 


1.69014638629 x 10" 15 


-3.4 x 10- 25 


19 





7.3 x 10" 14 


2.40241840942 x 10" 17 


-6.0 x 10~ 27 


21 





1.1 x 10" 15 


3.41486054743 x 10~ 19 


-1.0 x 10~ 28 


23 





1.7 x 10" 17 


4.85397236079 x 10~ 21 


-1.8 x 10~ 30 


25 





2.4 x 10~ 19 


6.89956364312 x 10~ 23 


-3.0 x 10~ 32 


27 





3.3 x 10~ 21 


9.8072207518 x 10" 25 


-4.9 x 10~ 34 


29 





4.2 x 10~ 23 


1.39402408398 x 10~ 26 


-8.1 x 10~ 36 


31 





5.2 x 10~ 25 


1.98150240103 x 10~ 28 


-1.3 x 10- 37 


33 





5.7 x 10~ 27 


2.81655949163 x 10~ 30 


-2.1 x 10~ 39 


35 





6.1 x 10~ 29 


4.00353154544 x 10~ 32 


-3.4 x 10" 41 


37 





6.2 x 10~ 31 


5.69072475939 x 10" 34 


-5.4 x 10" 43 


39 





5.9 x 10~ 33 


8.08894545219 x 10" 36 


-8.5 x 10" 45 


41 





5.0 x 10~ 35 


1.14978392551 x 10~ 37 


-1.3 x 10~ 46 


43 





3.6 x 10~ 37 


1.63433303287 x 10~ 39 


-2.1 x 10~ 48 


45 





1.7 x 10- 39 


2.32308384477 x 10~ 41 


-1.6 x 10~ 49 


47 











1.4 x 10" 50 



3.4 The exact solution of system (7) 

The approximate solution of system (7) is found. The aim of this section is to find an 
exact solution of (7), that is to say to find the function ip Q (x,t) in analytical form. For 
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arbitrary q G (0, 1) let us define the following sequence: 

n-l/2 

f = f { G IN : f 2n _, = T ^ rT ,/ 2n = }. 

The terms of the sequence f are proportional to Fourier coefficients of the Jacobi elliptic 
function cn (the elliptic cosine): 

7 00 / 7Z\ 27r 

cn(z, k) = - hn-i cos [i 2n ~ !) "J" J » where 7 = ^ , 2 G IR. (8) 

Let us clarify the notation and point out some properties of the elliptic cosine : 

• Basic periods of the doubly periodic function cn(^, k) are AK(k) and 2K(k) + 
2iK'(k), where K(k) is a full elliptic integral, K'(k) = K(k') and A;' = Vl — 



The parameter q in the Fourier expansion can be expressed in terms of elliptic 

— 7T — 

integrals: q = e a- . 



The Fourier-series expansion of the function cn(z, fc) does not include even har- 
monics. This expansion is valid in the following domain of the complex plane: 
—K' < 3m z < K', in particular, for z G IR. 

If z G IR and fc G (0, 1), then cn(z, k) G IR. 

The function cn(z, k) is a solution of the following differential equation: 

(2k 2 - 1) cn(z, k) - 2k 2 cn 3 (z, k). (9) 



d 2 cn(z,k) _ , , 



dz 2 

The latest property means that the infinite sequence of the Fourier coefficients for 
cn(z, k) is a solution of some infinite system of nonlinear algebraic equations. Let us find 
this system. On the one hand, it is clear from (8) that the Fourier-series expansion for 
the function cn 3 (z, k) is: 

cn 3 0, k) = ^ £ Fj 3) (f) cos y) ' where i = 1; 3 ' 5 ' • • • ' +°° ; 



f } (f) = 3// + 3/|/ 3i + } j 



( oo 2j-l \ 

6E(/ s 2 + / s /2 i+ J+3E/ s /2i- 
\ s+j s+j J 



+ 



oo oo oo oo j—2 J— 2 

~l~ ^ ^ ^ f s fpfj +s+p + 3 ^ ^ fsfpfs+p-j + E E fsfpfj-s 
s=l p=l s=l p=l s=l p=l 

s^j Pt^J s^j p^j 

(in all sums we summarize over only odd numbers). 
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On the other hand, from the differential equation (9) it follows that Fj 3 \f) is propor- 
tional to fj, with coefficients of proportionality depending on j: 

V, : ^ (fH pf^>4) /r (10) 

Thus, the sequence f is a nonzero solution of system (10) at all q G (0, 1). The following 
lemma proves the existence of the selected value of q. 

Lemma. There exists such value of parameter q G (0,1) that the 
sequence f becomes a real solution to system (7) with the 
real-valued parameter u 1 . 
Proof. Inserting the sequence f into system (7) : = fj and using system 
(10), we obtain: 



Vj : %(f) = ^ (3) (f)+/,(6E/n-32j% 

n=l 



f ~ „ 2 2(2k 2 - 1) , (\ 



' 3 \ Z-^/ J n 2 

n=l ' 



0. 



System (7) has a nonzero solution if and only if 

1 



1 256' 
< ~ (1 - 2fc 2 ) 

We have obtained the value of u^. The second equation of this system is equivalent 
to the following equation in parameter q: 

oo / n-1/2 \ 2 /i oo n \2 / oo n-1/2 \ 2 

3 S(rr^) -(3 + Stt^) +2 (£ttH = »■ (») 

This equation has the following solution on interval (0,1): 

q = 1.42142623201 x 10~ 2 ± 1 x 10~ 13 . 

Thus the lemma is proved. 

We can now construct the desired leading-order approximation for the function <p(x, t): 

ip (x, t) = A^cn(a(x — t),k) — cn(a(x + t), 

For arbitrary k £ (0, 1) this function is a real solution of equation (6). If a = ^f, then 
the periods of ip (x,t) in x and in t are equal to 2n. Using the Fourier- series expansion 
for the function cn(z,k) (formula (8)), we obtain the following expansion for function 

ip (x,t): 
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7 " 



<p (x, t) = 2Aj- / 2 n-i sm((2n - l)x) sin((2n - l)t). 



n=l 

If q = 1.42142623201 x 1(T 2 ± 1 x 10" 13 , then q is a solution of equation (11) and 
the sequence f is a real solution of system (7). The middle value of q corresponds to 
k = 0.451075598811, 7 = 3.78191440007 and a = 1.0576653982. All equations in system 
(7) are homogeneous, hence for these values of parameters, the sequence of the Fourier 
coefficients of the function tp (x, t) also is a solution of system (7), with 



2 

A 2 = 1.0983600974 A 2 . 



7 , 2 



1 64k 2 

Thus we have proved that the function 

(p (x, t) = A^cn(a(x — t),k) — cn(a(a; + 1), k)\, 

with k = 0.451075598811 and a = 1.0576653982 is such standing wave solution of equa- 
tion (5) that equation (6) has a periodic solutionis. 

3.5 The first approximation 

Now it is easy to find the periodic solution ip^x, t). Let us designate the Fourier coeffi- 
cients of the function <Pq(x, t) as D n f 



00 00 



(pl(x,i) = J2D nj sin(nx)sm(jt). 

n=l j=l 

Equation (6) gives the following result {b nn are arbitrary numbers): 

00 00 00 
i ) = Yl -2 I 3 2 sin ( nx ) sin(jt) + Kn sm(nx) sin(nt). 

n=l j=l J 71 n=l 

It should be noted that the function tp l (x, t) with arbitrary diagonal coefficients b nn is 
a solution of equation (6) and that all off-diagonal coefficients of <fx(x, t) are proportional 
to A 3 . 

3.6 The second approximation 

Let us consider the approximation to second order in e: 



g x 2 - g\ 2 = & 2 + {2UJ2 + Ul) dt 2 + 3V?l(X ' (12) 

If all diagonal coefficients of ifi(x, t) are zeros: Wn : b nn = 0, then Vj, n : bj n = — 6 nj - , 
and the function <fi(x, t)(fl(x, t) has no diagonal harmonics. Hence, selecting u 2 = — \uJ\ , 
we obtain a periodic solution to equation (12): 



x This result has been published in our paper pf 
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OO OO JJ oo 

<P2( X , f ] E E E - 2 ^ 2 sin ( nX ) Sin(jf) + k nn Sm(nx) SUl (fit), 
n=l j=l n n=l 



where 



if(a;,*) = 2cj x ^ h 3y> 1 (a;,t)¥>o( a; >*) = J2 J2 H nj sin(nx) sin(jt). 

°^ n=l j=l 

It should be noted that all diagonal coefficients h nn are arbitrary. 



4 Conclusions 

Using massless y? 4 theory as an example, we have shown that a uniform expansion of 
solutions for quasilinear Klein-Gordon equations can be constructed even in the principal 
resonance case. To construct the uniform expansion we have used the Poincare-Lindstedt 
method and the nontrivial zero approximation: the function 

tp (x, t) = A^cn(a(x — out), k) — cn(a(:r + out), fc)|, 

with k = 0.451075598811 and a = 1.0576653982. 

Thus, using the Jacobi elliptic function cn instead of the trigonometric function cos, 
we have put the principal resonance to zero and constructed, with accuracy 0(e 3 ), the 
doubly periodic solution in the standing wave form: 

(p(x, cut) = p (x, cut) + £(fi(x, cut) + e 2 Lp 2 (x, cut) + 0(e 3 ), 
with the frequency 

2 4 

cu = 1+^— A 2 e 1 -— A 4 e 2 +0(e 3 ) = l + 1.0983600974A 2 e-0.6031974518A 4 e 2 +C(£ 3 ). 

64P 8192A; 4 K J v J 

The question of constructing the uniform asymptotic series on the base of the obtained 
approximation remains open. The authors plan to investigate this question in future 
publications. 
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Appendix 



This program have been written in REDUCE 3.6. 

The program constructs system (7): Rjj = and finds a particular solution of this 
system. The function ip (x,t) has been denoted as phi(xx,tt). This function is a solution 
of equation (5). We assume that the function phi(xx,tt) is real and includes only odd 
harmonics: 

n 

phi(xx, tt) = a ( 2 J - !) sin((2j - l)xx) sin((2j - l)tt). 

3=1 

The number of unknowns (parameter n) can be selected arbitrary. For n — 10 this 
program gives result on computer with 16 Mbytes operating memory, for n = 20 this 
program gives result on computer with 128 Mbytes operating memory. 

Using the procedure "fourier" we expand equation (6) in the Fourier series and obtain 
the list of equations "listequa". The number of equations is 3n. The unknowns a (j) are 
a Fourier coefficients of some function, hence, we can assume that the sequence of the 
numbers c (j) = a (j)/a (1), satisfies the following condition: Vj > 1 : |c < 1. 

At the first step we assume that for j > 3, all c (j) = 0. From the first equation of 
"listequa" we find that C omega is a polynomial in c (3). We substitute this value of C omega 
in the second equation of "listequa". As result we obtain that this equation is a cubic 
equation in c (3). We solve this equation, using the standard procedure SOLVE. One of 
its solution must be real and we select this solution as value of c (3). 

At the second step we assume that c (5) is arbitrary { clear c (5) }. From the first 
equation of "listequa" we obtain that C omega is a polynomial in c(5), substitute this value 
of C omega in the third equation of "listequa" and solve it, determining c (5). 

At the third step we substitute new values of c (5) and C omega in the second equation 
of " listequa" , so we calculate the variable "equa". If "equa" is more than 5 = 10 -11 , then 
we resolve the second equation of "listequa" and repeat the second step. Alternatively, if 
the obtained value of "equa" is less than or equal to 5 = 10 -11 , then we assume that the 
three leading equations of "listequa" are solved and seek c(7) from the fourth equation of 
"listequa" and so on. Once the coefficient c (2n — 1) has been found the program finishes 
to work. 



out "rjj .res" ; 
in fourier$ 
n:=10$ 

operator a,c,phi$ 
depend phi,xx,tt$ 

% Expand the unknown function phi(xx,tt), which is a solution of 
% equation (5), and equation (6) in the Fourier series. 

phi(xx,tt) :=for j:=l:n sum a(2j-l)*sin((2j-l)*xx)*sin((2*j-l)*tt)$ 
equationl :=2*wl*df (phi(xx,tt) ,tt ,2)+ 
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f our ier (four ier (phi (xx , tt ) **3 , xx) , tt ) $ 

listequa: =0$ 

°/ The frequency w_l is proportional to a(l)~2. 
wl :=C_omega*a(l)**2$ 

for j:=2:n do a(2*j-l) :=c(2*j-l)*a(l)$ 
equationl :=equationl$ 

% Construct the system (7): R_{j j}(c)=0. 

for k:=l:3*n do 

li st equa:=append (listequa, 

{16*df (equationl , sin( (2*k-l) *xx) , 

sin((2*k-l)*tt))/a(l)**3}) ; 

°/ Find a real solution of this system. 

on rounded$ 
C_omega:=0$ 

°/ C_omega must be connected variable, because I want to clear it 
% without message : "WARNING . . . " . 

for k:=2:n do 

« for j:=k:n do c(2*j-l):=0; 
clear C_omega, c(2*k-l); 
f irstequa:=f irst (listequa) ; 

C_omega : =C_omega-f irstequa/ df (f irstequa , C_omega) ; 
equa: =part (listequa, k) ; 

°/ All equations, which we solve using the procedure SOLVE, are 

°/„ cubic equations, hence, they have at least one real root. We assume 

°/ that the absolute value of one of real roots is less than unit. 

°/„ We solve the equation and select this real solution 

% as value of c(2*k-l) . 

solve_equa: =solve (equa, c (2*k-l) ) ; 
while solve_equa neq {} do 
« c(2*k-l) :=part(solve_equa, 1,2) ; 

if( c(2*k-l)=sub(i=-i,c(2*k-l)) and abs(c(2*k-l))<l) 
then solve_equa:={} 
else solve_equa: =rest (solve_equa) ; 

»; 

test :=0; 
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while (test=0) do 
« test:=l; 

clear C.omega; 

f irstequa: =f irst (listequa) ; 

C.omega : =C_omega-f irstequa/df (f irstequa , C_omega) ; 
for j :=2 :k do 

« 

equa:=part(listequa, j) ; 
if (abs(equa)>10**(-ll)) then 
« clear c(2*j-l) ; 
test :=0; 

equa :=part (listequa, j) ; 
solve_equa : =solve (equa, c (2* j -1) ) ; 
while solve_equa neq {} do 

« c(2*j-l) :=part(solve_equa,l,2) ; 

if (c(2*j-l)=sub(i=-i,c(2*j-l)) and abs(c(2*j-l))<l) 
then solve_equa: ={} 
else solve_equa: =rest (solve_equa) ; 

»; 

»; 
»; 
»; 
»; 

°/ Write the obtained result . 
c(l) :=1$ 

for j:=l:n do write "c(" ,2*j-l ,"):=" , c(2*j-l); 

for j:=2:n do write "c(" ,2*j-3, ")/c(" ,2*j-l, ") := ", c(2*j-3)/c(2*j-l) ; 
write "C.omega: =" , C_omega; 

for j :=l:3*n do write "R(" ,j ,j ,"):=" , part (listequa, j ) ; 
quit ; 

The procedure "fourier" expands polynomials of sin(a:) and cos(x) in Fourier series. 

procedure f ourier (FF,X) ; 

% This procedure constructs Fourier-series expansions for polynomials 
% of sin(x) and cos(x). 

begin 

scalar F; 

for all a, b such that df(a,x) neq and df(b,x) neq 
let 
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cos (a) *cos (b) = (cos (a-b) +cos (a+b) ) /2 , 
sin(a)*sin(b)=(cos(a-b)-cos(a+b) )/2, 
sin(a)*cos(b)=(sin(a-b)+sin(a+b) )/2, 
sin(a)**2=(l-cos(2*a))/2, 
cos(a)**2=(l+cos(2*a))/2; 

F : =FF ; 

for all a, b such that df(a,x) neq and df(b,x) neq 
clear 

cos (a) *cos (b) , 
sin(a) *sin(b) , 
sin(a) *cos (b) , 
sin(a)**2, 
cos(a)**2; 

return F; 

end; 
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